Superharmonic long-range triplet current in a diffusive Josephson junction 
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We study the Josephson current through a long ferromagnetic bilayer in the diffusive regime. 
For non-collinear magnetizations, we find that the current-phase relation is dominated by its second 
harmonic, which corresponds to the long-range coherent propagation of two triplet pairs of electrons. 
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The interplay between ferromagnetism and supercon- 
ductivity in hybrid structures is an active area of research 
0, In a homogeneous ferromagnet (F) adjacent to a 
conventional superconductor (S), superconducting corre- 
lations between pairs of electrons are induced both in the 
spin singlet channel and the triplet channel without spin 
projection along the magnetization axis. This proxim- 
ity effect is short-ranged due to the dephasing of elec- 
trons with opposite spins induced by the ferromagnetic 
exchange field. On the other hand, in the presence of a 
non-collinear magnetic configuration, a long-range prox- 
imity effect can be induced in the triplet channels with 
parallel electron spins, as no such dephasing occurs in 
that case 0, 01 

Early proposals to observe this effect suggested mea- 
suring the critical current in a Josephson junction 
through a ferromagnetic trilayer with non-collinear mag- 
netizations [H, Q . In this geometry, the external F lay- 
ers convert singlet Cooper pairs into triplet pairs which 
have a non- vanishing projection onto the channels with 
parallel electron spins along the (tilted) magnetization 
of the central F layer, and thus may propagate coher- 
ently over long distances. A maximal critical current 
with amplitude comparable to that of a normal metal- 
lic Josephson junction with the same length is obtained 
when the external layers have a thickness comparable to 
the ferromagnetic coherence length, £p, and the mag- 
netizations in successive layers are orthogonal. (Simi- 
larly, long-range triplet correlations can be induced by 
spin-active interfaces between a ferromagnet and the su- 
perconducting leads 04^-) Indeed, recent experiments 
observed a strong enhancement of the Josephson current 
through a ferromagnetic multilayer when the layers were 
ordered non-collinearly [l"ol - fl2j . 

Here we arc interested in the question whether three 
layers (or, equivalcntly, spin-active interfaces on both 
sides of the junction) are necessary to observe a long- 
range triplet Josephson current. Recently it has been 
shown theoretically that a long-range triplet proximity 
effect may also develop in ballistic bilayer ferromagnetic 
Josephson junctions with non-collinear magnetizations 
fl3j . In this case, a superharmonic Josephson relation 
is generated by the long-range propagation of an even 
number of triplet pairs which may then recombine into 




FIG. 1: Setup of the junction. The superconducting leads 
are coupled through two ferromagnetic layers in series, with 
thicknesses cIl and oIr, respectively. The magnetizations of 
the layers are tilted with respect to each other by an angle 9. 



singlet Cooper pairs. It is important to know whether 
this effect is robust to disorder which is believed to be 
present in the experiments, i.e., whether it exists in dif- 
fusive systems as well. In this work, we show that this 
is indeed the case. We find that the amplitude of the 
critical current determined by the second harmonic of 
the Josephson relation decays algebraically with length. 
The maximal current is smaller than in trilayers due to 
fact that the triplet pairs have to recombine into singlet 
Cooper pairs, a process that takes place only within a 
distance ~ £f near the F/S interface. 

We consider a Josephson junction formed of two super- 
conducting leads contacted through two ferromagnetic 
layers in series, with thicknesses and d^, respectively, 
cf. Fig. [T] The F layers are assumed to have the same 
properties, but different orientations of their magnetiza- 
tions. As we are interested in the long-range proximity 
effect, moreover, we assume that the right layer of length 
dfi is much thicker than the diffusive ferromagnetic co- 
herence length, £f = y/D/h, where D is the diffusion 
constant and h is the amplitude of the exchange field. 

Within the quasiclassical diffusive theory [14| , the equi- 
librium supercurrent flowing through the junction can 
be expressed through the quasiclassical Green function 
g, which is a 4 x 4 matrix in the particle/hole and spin 
spaces, and obeys normalization conditions g 2 = 1 and 
Tr g — 0. It solves the nonlinear Usadel equation 



- Dd v {gd y g) + [(w + ih(y).a)T z ,g] = 0, (1) 
where y is a coordinate along the junction, Tj and <jj 



2 



= x, y, z) are Pauli matrices in the particle/hole and 
spin spaces, respectively, and oj = (2n + 1)ttT (n integer) 
is a Matsubara frequency at temperature T. The orien- 
tations of the exchange fields are characterized by a tilt 
angle, 



h(y) = h[s'm 6(y) x + cos 9(y) z] 



(2) 



where 9{y) = 9 for — di < y < and 9{y) = for 
< y < d,R. Note that the orbital effect is neglected in 
Eq. (fl]) by assuming that the magnetic flux through the 
junction is much smaller than the flux quantum. 

Equation (JTJ) has to be supplemented by boundary con- 
ditions at each interface. We assume that the interface 
resistance between the two F layers is much smaller than 
the resistance of each layer. Therefore, g and its deriva- 
tive are continuous at y = 0. Moreover, assuming a rigid 
boundary condition at the interfaces with the leads as 
well as negligible resistances of the F/S interfaces, we im- 
pose the continuity conditions g(— d^) = g^ and g{dp) = 
g R , where g k = [ujt z + A(cos (p k T x - sin tp k T y )}/\/u} 2 + A 2 
(k — L, R) are the Green functions of the leads. Here A 
is the amplitude of the superconducting gap and ip^ its 
phase. For convenience, we choose tfL = if and ipn = 
such that tp denotes the superconducting phase differ- 
ence. 

The supercurrent is then related to the Green function 
g through 



I 



ivDST ^2 lmTr b~ 3 gd y g], 



(3) 



uj>0 



where v is the density of states (per spin) at the Fermi 
level and S is the cross section of the junction. Due to 
current conservation, Eq. §3§ may be evaluated at any 
position along the junction. 

To proceed further, we introduce the parametrization 



(4) 



which automatically satisfies the normalization condi- 
tions. Here, the anomalous functions F and F = a y F*a y 
are 2x2 matrices in spin space and odd in A. With the 
parametrization the Usadel equation ([1]) takes the 
form 




Dd„ 



Vl-FFdyF - F d y Vl-FF 

+ {w + ih(y).cr,F} = 0, (5) 

while the current ([3]) simplifies to 

I = 2-izevDST ^ 1mTr[Fd v F]. (6) 



uj>0 



Note that, as F and F arc odd functions of A, the current 
is even in A. 



At temperatures slightly below the superconducting 
critical temperature T c , the gap vanishes as A cx 
\/T c {T c — T). Thus, one may solve Eq. §5§ perturba- 
tively around the normal state solution, F = 0. To this 
end, we expand F in the small parameter (l5| A/w, i.e., 
F = (A/o;)FW + (A/cofF^ + and solve Eq. © 
order by order. To obtain the current up to the fourth 
order in A/T c , it is sufficient to compute F^ and F^ 3 \ 

Let us start with the leading order. Then, Eq. (|5j) 
yields the linear differential equation, 



D8 2 y F^ 



2loF {1) +i{h(y).(T,F^} =0. (7) 



Upon performing a unitary transformation, F^'(y) = 
e ia v 0( y )/2 F (i)^ e -i<7 y e(y)/2^ thc gcnera i solution in each 

layer takes the form F^ = Fq Pcrx+Fz 1 cr z , where 



with s 



x-(i) = a l/r - 

L/R,s - ^-s e 

±,x and F± = 



P=V _|_ gL/R e p a y 



(8) 

Here F { ± ] 



with p± = y/2(u> ± ih)/D correspond to the short-range 
singlet and triplet correlations, while F^ with p x = 
y/2ui/D corresponds to an equal superposition of the 
long-range triplet correlations. For typical ferromagnets, 
h > T c , therefore, p± « (1 ± 

The coefficients A^ R and Bs that enter Eqs. ((HJ 
are determined by the boundary conditions. Assuming 



dn > £f, we find 
AR e -P- dR = ll'j < <■ - ! - : = ie lv a.. 
small corrections in c?r/£f) with 



- p±dR , and 



(up to exponentially 



a+ = a_ = 



a- 13 



a 



tan 9 
2 smhp x d 



t, (9a) 
Im [(p+d L + l)o+] , (9b) 



and 



o 



cos 2 -e p ~ dL 



sm 



2 cos 9 



(p-di 



coshp_d L smhp x d R 

1) r-r — ~ , (10a) 

Sinn a 



sin — (coshp + (ii — i sinhp+di) 

sin 2 6> coshp + di sinhp^dfl 

-7; n(P- d L + V r-r ~ , 10b) 

2cos# smhp x a 

where d = d^ +d R is the total length of the junction. We 
do not give the explicit results in the left layer because, 
in the following, we choose to evaluate the current (JB]) in 
the right layer. 

Note that the ratio d^/^F may be arbitrary, while 
Eqs. ((9]) and (fT0|) have been simplified with the assump- 
tion di < £jv, where £w = sj D/ (2nT c ) (£ N > £ F ) is the 
normal coherence length close to T c . For a non-collinear 
configuration of the layers, long-range correlations are 
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FIG. 2: Dependence of the short-range first harmonic h 
of the current-phase relation on the thickness gIl/^f and the 
angle 6 at fixed dn/^F = 5. Here h is measured in units of 
ho = 7vGA 2 d/(V2eT c ^ F )e- d/ ^. 



present (F^ ^ 0). Nevertheless, they do not contribute 
to the first harmonic of the current-phase relation, as we 
show now. 

Inserting the solution (jHJ for F R into Eq. ([6]). we ob- 
tain the first harmonic I\ of the current-phase relation, 



= /isintys, where I\ oc Re[p+e 
the form 



-p+d R 



a+\ 



ttGA 2 d 
h = —f= 

V2eT c £f 



d/£ FT I ^R_ 
1 Cf'Cf' 



It takes 



(11) 



where G = 2e 2 vDS/d is the conductance of the junction 
in the normal state and Ti is a scaling function. 

The exponential reduction factor e~ d ^ F in Eq. (fTTj) 
shows that I\ is always short-ranged. The scaling func- 
tion Ti is shown in Fig. [2] as a function of the length 
of the left ferromagnet and the tilt angle between the 
magnetizations. Its maximal absolute value is of the or- 



der of 1 and its sign oscillates, thereby displaying transi- 
tions between 0-states (when 1\ > 0) and 7r-states (when 
h < 0). Such oscillations have indeed been observed in 
S/F/S junctions 

Simple analytic expressions for the function Ti in 
Eq. can be obtained in various regimes. For a par- 
allel alignment, = 0, the results for a long (d ^> £f) 
monodomain S/F/S junction are reproduced with Ti ~ 
sin(d/^F + tt/4) [l9(. Moreover, as long as g?l -C £f, this 
result remains valid for arbitrary angles. In the opposite 
limit, cLl ^S> £,f, wc find Ti ~ y/2 cos(dL / £f) cos(cZr/£_f) 
for 89 <9 <n-89 and Ti ~ V2cos[(d L -d R )/( F } for an 
antiparallel alignment, 9 = ir. Here 89 ~ ^/min(dL) dn) 
is the crossover scale over which Ti varies smoothly be- 
tween its respective values at 9 = 0, at intermediate an- 
gles, and 7T. These results generalize those of Refs. 20l 21\ 



obtained at dr = dn 22 



We now turn to higher-order contributions in A/T c to 
the current-phase relation. The fourth-order terms con- 
sist in a small (short-range) correction to the amplitude 
1 1 of the first harmonic, as well as a second harmonic 
oc sin(2y?). Here we concentrate on the second harmonic. 
In collinear structures, its amplitude li is short-ranged 
with a suppression factor e~ 2d ^ F . Thus, it becomes im- 
portant only near 0/-7r-transitions [23l. l24j when I\ van- 
ishes. In non-collinear structures, however, I2 is long- 
ranged as we show now. 

For concreteness, we will assume dn 3> c?l and 
evaluate the current ((6J in the right layer. Long- 
range contributions originate from the propagation of 
parallel-spin triplet pairs described by F x , namely 1% oc 

lmTr[F^d y F R % + F^d y F^ x ). Using Eq. ©, we find 

(3) 

that F R obeys the equation 



J 



2^ x = ^d y {[(F^r + F^l 



F (l) Q pW 

r R, + u V r R,+ 



d f (1) 

U V C R- 



7 (1)1 

R,x J 



(12) 



r 



Being interested in the second harmonic, we only need 



(3) 

to retain contributions to F R x that are proportional 
e~ ltp . Thus, we may decompose F^ x = f x 3) e~ tv + $ 

(3) 

where Q R x does not contain any contributions oc 

(3) 

The function fx then obeys a differential equation ob- 
tained from Eq. (fT2j) by keeping on the r.h.s. the terms 
oc e~ lip only. In particular, such terms arise from the 



to 

(3) 

R,xi 

-itp 



contribution F^l^F^^dyF^ 1 ^ yielding 

, 3)^ 9 ,„ f (3)__^ 3 
J yJx 



^/j 3 > + 2^> = -^^a. 



where we used that £f *C £,n- A general solution that 



satisfies Eq. (fTBf together with the boundary condition 
fx [d R ) = at the right F/S interface reads 



,(3) = P^f A 3 



(14) 



2 u 
+Csinhp x (d R - y) 



Here C is a constant which should be determined from 
the boundary condition at the interface between the left 
and right ferromagnet. It turns out, however, that it does 
not contribute to the current and, thus, will not need to 
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be determined. We obtain 



u»0 



I 2 =^eiySDTj2^P 2 AFal 
which evaluates to 



(15) 



!92eTid R 2 Uf' 



384d| 2d R /t. 



(16) 

The second harmonic is long-ranged: it depends on 
cLr/S^f as a power law. The suppression factor £f/^r 
is due to the conversion of parallel-spin triplet pairs into 
singlet Cooper pairs which takes place on a distance £p 
from the F/S interface only. The dependence of the func- 
tion T2 that appears in Eq. (|16[) on the tilt angle and 
the thickness of the short layer is shown in Fig. [3] T2 
is maximal for lengths ~ £p and angles 9 ~ ir/2. As 
expected, it vanishes for collinear structures, 9 — > 0,7r of 
arbitrary length. In non-collincar structures, it simplifies 
to T 2 ~ (<W£f) 4 sin 2 9 for a short F layer, d L -C £f, 
while it vanishes exponentially 



4 sin 



(cosf^ 



cos 



2 9 



sin 



(17) 



for a long F layer, dL S> £f- Here ag is only relevant in 
a vicinity (5$ near 9 = 0,7r, where ao = 2 and a T = 1, 
respectively. 

As in non-collinear structures, the second harmonic I2 
is long-ranged (in contrast to the first harmonic 2i), we 
expect it to dominate as soon as d R exceeds a few times 
£p. In particular, for an optimal thickness d^ ~ we 
find the ratio I 2 /h ~ {A/T c ) 2 ^ F /d R ) 2 e dR /^ close to 
T c . The effect is expected to be robust at lower temper- 
atures, further enhancing the ratio I2/I1 when A ~ T c , 
More generally, odd and even harmonics will be short- 
ranged and long-ranged, respectively. However, as triplet 
pairs need to recombine into singlet Cooper pairs on a 
distance £p from the F/S interface, we conjecture that 
the amplitude l2n of the even harmonics oc $,m.(2rnp) will 
contain a small factor (^/(ifl)", sec Eq. (fT6")) for the case 
n = 1. Thus, the current-phase relation could be domi- 
nated by the second harmonic at all temperatures. 

Let us now discuss the measurability of our predic- 
tion. The long-range proximity effect has been observed 



in trilayers Hj, 11 



where already the first harmonic is 
long-ranged. The amplitude of the second harmonic for 
bilayers predicted in our work is smaller only by a fac- 
tor ^f/<1r as compared to the first harmonic in trilay- 
ers. Thus we believe that it should be well within the 
sensitivity of present-day experiments. Its specific phase 
dependence may be detected by a direct measurement of 
the current-phase relation [iH or through the appearance 
of fractional Shapiro steps in the current-voltage charac- 
teristics under microwave irradiation 12611. 




FIG. 3: Dependence of the long-range second harmonic 
I2 of the current-phase relation on the thickness dhj^F 
and the angle 8. Here I2 is measured in units of I20 = 
7rGA 4 £ F /(192eT c 3 dK). The plots b) and c) are taken along 
the cuts indicated by the blue lines in a), at 9 — 0.597T and 
di/^F = 1.26, respectively. Both cuts include the maximal 



In conclusion, we predicted that the current-phase 
Josephson relation through a long diffusive ferromag- 
netic bilayer with non-collinear magnetizations is dom- 
inated by a superharmonic contribution ex sin(2tp). The 
second harmonic can be viewed as the minimal Joseph- 
son current 2JJ that can flow between the conventional 



even-frequency superconductor in one lead and the effec- 
tively odd-frequency superconductor 0, [28| generated by 
the long-range triplet proximity effect at the extremity 
of the ferromagnetic bilayer attached to the other lead. 
Measuring the dependence of the Josephson current on 
the thicknesses of the ferromagnetic layers and the angle 
between their magnetizations would provide further evi- 
dence for the long-range triplet proximity effect. Similar 
experimental studies have been performed for the criti- 
cal temperature of a thin superconducting film in con- 
tact with a ferromagnetic spin valve (29| . Furthermore, 
detecting the 7r-periodicity of the current-phase relation 
through phase-sensitive measurements would be a strong 
indication of the odd-frequency nature of the long-range 
proximity effect. 
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